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ABSTRACT 

The strain induced by lattice mismatch at the interface is responsible for the different 
value of the band discontinuities observed recently for the AlN/GaN (AIN on GaN) and the 
GaN/AlN (GaN on AIN) polar (0001) interface. We present a first-principles calculation of 
valence band offsets, interface dipoles, strain-induced piezoelectric fields, relaxed geometric 
structure, and formation energies. Our results confirm the existence of a large forward- 
backward asymmetry for this interface. 

INTRODUCTION 

A reliable determination of the valence-band offset (VBO) at the (0001) polar interface 
between wurtzite AIN and GaN is still missing. The few experimental investigations avail- 
able JI], 3 are in mutual disagreement, and theoretical studies refer either to zincblende 
interfaces [^, or artificially lattice-matched wurtzite interfaces [Q. The latter approxima- 
tion leads to a less accurate determination for the VBO, and cannot pick up the possible 
forward-backward asymmetry characteristic of lattice-mismatched interfaces. In the case 
of the AlN/GaN interface, lattice mismatch amounts to 2.5 %, and may cause a very large 
asymmetry. This asymmetry has not yet been clearly determined experimentally (it was 
not even found in early experimental work 0), being hidden by the large uncertainties in 
the measured data. Even the best experimental investigations available face two kinds of 
problem: (i) the determination of the core level alignment with the valence-band maximum 
(VBM) is obtained indirectly using theoretical estimates of the VBM position, which (as 
underlined by Vogel et al. p) is affected by large systematic errors; (ii) the existence 
of strong polarization fields in both the substrate and the overlayer tends to modify the 
apparent value of the VBO deduced from the core-level shift measurements. 

The present ab-initio investigation includes all strain and relaxation-induced effects, 
and overcomes the difficulty in VBO determination due to polarization fields by the use 
of a novel charge-decomposition technique. An estimate of the formation energy of the 
interfaces studied is also given. 

BULK PROPERTIES 

Valence-band offset calculations at lattice-mismatched interfaces require the evaluation 
of the band structure energies for the bulk crystals in equilibrium, and subjected to biax- 
ial strain. The calculations are done using density-functional theory in the local-density 
approximation (LDA) to describe the exchange-correlation energy, and ultrasoft pseudopo- 
tentials for the electron-ion interaction. Plane-wave basis sets up to 25 Ry, and 24 
special k-points are found to give fully converged values for the bulk properties. Since the 
properties of GaN are affected by Ga 3d states our Ga pseudopotential includes 3d 
electrons in the valence. This yields very good structural bulk parameters (see below). 



Table I: Predicted structural parameters and valence band maxima for equilibrium and 
strained AIN and GaN. 



Material 


AIM 


AIM 


P qM 


P qM 


Substrate 




GaN 




AIN 


a 


5.814 


6.04 


6.04 


5.814 


c/a 


1.619 


1.51 


1.6336 


1.73 


u 


0.38 


0.3927 


0.3761 


0.3653 


-E^strain (eV) 




+0.179 




+0.155 


^VBM (eV) 


-0.16 


0.09 


-4.90 


-4.69 



In wurtzite crystals, the determination of the atomic structure at a given lattice constant 
a implies the calculation of the c/a and u parameters. The equilibrium c as been determined 
by fitting with a polynomial the total energy computed for six different values of c, with u 
being determined for each value of c/a via minimization of the Hellmann-Feynman forces, 
with a threshold of 10~^ Hartree/bohr. The calculated structural parameters are given 
in Table I. The structural parameters of AIN and GaN behave similarly under strain 
[A(c/a)/(c/a) ~ 7%, Au/u ~ 2%] with similar total energy variations. Instead, the effect 
of strain on the valence-band edge is very different. A rationale for this difference is that 
the AIN (GaN) band edge is a singlet (doublet) formed by the hybridization along the c-axis 
(in the a-plane) of N 2s orbitals with Al (Ga p^y) states, so that biaxial compression 
pushes the edge upward in GaN and downward in AIN. 

BAND OFFSET 

As pointed out by Baldereschi et al. 0, the valence-band offset AE„ may be split in 
two terms: 

AE, = AEvBM + AVei. 

The first contribution Aii^vBM is the difference between the valence-band edge energy in 
the two bulk materials, each edge being referred to the average bulk electrostatic potential. 
The second contribution, the potential lineup AVei, is the drop of the macroscopic average 
of the electrostatic potential across the interface. The latter term requires a selfconsistent 
calculation of the electronic density distribution for the real interface system. Our interface 
has been modeled using a (GaN)4/(AlN)4(0001) superlattice (see below), both ideal and 
fully relaxed. The material being grown epitaxially on the chosen substrate, has been 
pre-strained to have the same a lattice constant as the substrate. 

The lineup term is customarily obtained by solving the Poisson equation for the macro- 
scopic average of the charge density, neutralized by a suitable distribution of gaussian 
charges centered on the ion sites. The potential drop across the interface is usually cal- 
culated as the difference of potential values in bulk-like regions inside the two interfaced 
materials. This turns out to be non-trivial for a system such as the present one, in which 
the existence of polarization fields in the equilibrium bulk makes it impossible to define 
asymptotic values for the electrostatic potentials. The existence of such fields, moreover, 
limits the maximum length of our slab. Indeed, beyond a certain critical thickness the drop 
of the potential inside each slab would make the system metallic, with a related transfer 
of charge, which would spoil the exact determination of the lineup term. Our choice of a 



16-atom supercell is a compromise between the need to have an insulating system, and at 
the same time to avoid a spurious couphng of the interfaces at the sides of the slabs. Tests 
were performed in supercells of up to 40 atoms. In Fig. 1, we show the macroscopic average 
of the charge density and of the electrostatic potential. The potential drop is inextricably 
linked to the polarization fields within the AIN and GaN bulks. We have circumvented 
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Figure 1: Supercell electron density and electrostatic potential. Electron density has been 
compensated in the two bulks by a distribution of gaussians placed at the atomic sites. 

this problem by employing a new method. The basic idea is that at the polar AlN/GaN 
interfaces, the existence of polarization fields reveals itself by an accumulation of charge in 
the form of a monopole distribution whose density is proportional to the difference between 
the polarizations inside the two interfaced bulks. On top of this monopole term, we have 
the traditional dipole term representing the local charge transfer across the interface. This 
dipole term is the quantity we are interested in, as the band offset is by definition related to 
the dipolar part of the potential drop. Since the monopole contributions are related to the 
polarization fields, they must be equal and opposite for the two (geometrically inequivalent) 
interfaces in our AlN/GaN superlattice. To filter out the monopole term we superimpose 
the two interface distributions by folding them around a plane placed halfway between the 
two junctions. We define the dipole term n*^*^ as the average of the superimposed charges, 

n'^'Piz -Zo) = ^ [n{z - Zo) + n{zo - z)] , 

where 2; is a coordinate along the c-axis, zq the plane position and n the macroscopic average 
for the charge density. The monopole term fi^°^° is just the difference between the dipole 
term and the total macroscopic charge: 

n™{z) =n{z)-n'^'^{z). 



Such a decomposition allows a determination of the polarization charges and dipole terms 
which is nearly independent of the position of the folding plane. Fig. 2 shows the decompo- 
sition for the AlN/GaN interface. The decomposition reveals the origin of the asymmetry 
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Figure 2: Decomposition of the macroscopic average of the electronic density (dotted line) 
into monopole (solid) and dipole (dashed) terms. Such a decomposition allows the deter- 
mination of the lineup term AVez from the solution of the Poisson equation (dot-dashed) 
of the dipole term. 

in the total charge distribution, and at the same time it enables us to evaluate the lineup 
term. In Table II we report the values for the VBO obtained via this decomposition. There 



Table II: Valence-band offset Ai?„, potential lineup AVei, relaxation energies Erei, monopole 
charge densities aint and electric fields E in the ideal and relaxed the AlN/GaN (0001) 
interface. 



Interface 


AlN/GaN 


GaN/AlN 


units 


structure 


ideal 


relaxed 


ideal 


relaxed 




AE, 


0.29 


0.20 


1.00 


0.85 


eV 


AV,i 


5.28 


5.18 


5.52 


5.36 


eV 


<^int 


0.029 


0.014 


0.022 


0.011 


C/m2 


E 


32.7 


15.6 


24.4 


12.9 


10^ V/m 



is a very large forward-backward asymmetry of 0.65 eV between AlN/GaN and GaN/AlN 
interface VBOs. This is only marginally due to the lineup term (contributing 0.18 eV), its 
main component being the band structure term (0.47 eV). The relaxation is responsible for 



comparatively small deviations of ~ 0.1 eV from the ideal-interface values. The relaxation 
pattern is characterized in both cases by a contraction of the Al-N axial interface bond (~ 
-0.04 a.u.) and an expansion of the axial Ga-N bond (~ +0.02 a.u.). 

POLARIZATION 

Supercell calculations are not the only way to obtain the interface charge density aint- 
As shown in Ref . [|l^ , given the polarization Pi and P2 and the dielectric constants 61 and 
62 of the component materials, aint is given by 



a,„i = ±2(P2-Pi)/(£i + £2) 



(1) 



We have calculated the macroscopic polarization for 



in periodic boundary conditions. 

equilibrium and strained GaN and AIN via the Berry phase technique of Ref. 



The 



(high-frequency) dielectric constants of AIN and GaN have been calculated using the rela- 
tion 

APt = £ooAPl, 

where APp is the (so-called transverse) polarization change induced by a small cation 
sublattice displacement in the bulk in zero field, and APl is the (so-called longitudinal) 
polarization change due to a uniform displacement of few cation planes in a periodic bulk 
supercell. In the latter, a depolarizing field is present due to the periodic boundary con- 
ditions. As a by-product of our calculations, we obtained the Born effective charges for 
AIN and GaN which, as expected for highly polar semiconductors, are quite close to the 
nominal ionicity. The results are shown in Table III. Substituting P^^ and e in Eq.[^ we 

Table III: Polarization in GaN and AIN: electronic Pei and total Ptot polarization, derivative 
of the latter with respect to u, Born effective charge, and dielectric constant are shown. 



System 


a 


Pel 


Ptot 


dP/du 


Z* 


^00 




bohr 


C/m^ 


C/m2 


C/m" 


e 




AIN 


5.814 


-0.178 


-0.0812 


10.51 


2.69 


4.59 


AIN 


6.04 


-0.492 


-0.1712 


9.95 


2.74 


4.64 


GaN 


5.814 


+0.223 


+0.0343 


10.69 


2.74 


5.27 


GaN 


6.04 


-0.0511 


-0.0308 


9.88 


2.72 


5.52 



obtain for the AlN/GaN (GaN/AlN) interface a monopole density of 0.028 (0.023) C/m^ 
in embarrassing agreement with the outcomes of the supercell calculations. This proves 
directly that the sources of the internal fields in the interface system are the charges accu- 
mulated at the interface by the bulk polarization effects. Also, the latter finding provides 
an a posteriori justification of our somewhat ad-hoc charge decomposition procedure. 

FORMATION ENERGIES 

An important issue for the present selfconsistent calculation is the evaluation of the 
formation energy for the AlN/GaN interfaces. Contrary to the case of non-polar interfaces, 
for the wurtzite (0001) system it is impossible to build a superlattice with symmetric 
interfaces. This means that only the average value of the formation energy for the two 



interfaces can be obtained from a total-energy calculation. We define the average formation 
energy for the AlN/GaN interface as 



^ 2A 



E 



tot 



^Ga^GaN 



where /j,^ are the total energies per pair of GaN and AIN, the number of Ga and Al 
atoms, the supercell total energy and A its cross-sectional area. A reliable determi- 
nation of Ef requires equivalent k-point sampling for bulk and interface calculations. This 
is easily accomplished if the interface is lattice-matched. In the present case, the supercell 
length is not simply an integer multiple of the bulk unit cell of either constituent material. 
This means that an exact equivalence between k-point meshes cannot be achieved. A good 
approximation for Ef can however be obtained by defining, for each component material, 
an auxiliary bulk cell having the same lattice constant a, and an axial length c being a 
sub-multiple of the supercell length /. This value in the present case is just the average of 
Cain and ccaN- It is then possible to downfold exactly the supercell mesh into the auxiliary 
bulk cell. The next step is to uniformly scale the k-points coordinates to adapt the mesh 
to the real value of c. We should point out that the accuracy of this procedure (compared 
with an exact computation of the energy integral over the IBZ) increases with the number 
of points in the mesh. It is therefore possible to find a suitable mesh to accomplish any 
required accuracy. The results for the formation energies reported in Table |I^ have been 

Table IV: Average formation energy for the AlN/GaN (0001) interfaces. 



Interface 


AlN/GaN 


GaN/AlN 


units 




ideal 


relaxed 


ideal 


relaxed 




Ef 


-3.4 


-16.4 


+11.7 


-6.2 


meV 



obtained using a 6-point Chadi-Cohen mesh ||T2[ in the supercell, which when downfolded 
in the auxiliary cell produces 24 special points. We estimate the k-point sampling error 
in the formation energies to be ~ 10 meV. It should be noted that such low formation 
energies are not surprising when compared with the results obtained by Chetty et al. |]r3 
for GaAs/AlAs(lll) interfaces. 
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